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Abstract 

The paper proposes a different method of solving a simplified version of the Black-Scholes 
equation. In the first part of the paper, the Black-Scholes equation is transformed into ordinary 
differential equation to get a solution similar to the solution of the Euler equation. The second 
part of the paper focuses on partial differential equation. Separation of variables method is used 
to solve the Black-Scholes equation. Plots corresponding to put and call options are also given. 


1 Introduction 

Differential equations have a great variety of applications in different fields of science such as 
engineering, physics, biology, pharmacokinetics (Li et al. (2014)). Yet, there are only a few of their 
applications in economics or finance. Particularly, well-known models involving differential equations 
are only economic growth model and Black-Scholes equation. The latter one will be discussed in 
the paper. In 1977, Myron Scholes together with Fischer Black were awarded a Nobel Prize in 
economics for the formulation of stock options formula through “new method of determining the value 
of derivative” (Jarrow (1999)). So, Black-Scholes model deals with one of the most important issues 
in quantitative finance pricing of options (Rodrigoa and Mamon (2006)). This model has significant 
implications both theoretical and practical since finance plays a great role in economies around the 
world (Bohner and Zheng (2009)). 

2 Ordinary Differential Equation 

2.1 Background information and underlying assumptions 

In practice, Black-Scholes model of option pricing was applied to various “commodities and payoff 
structures” (Jodar et al. (2005)). Black-Scholes model is widely used for American options as well as 
for European options. Therefore, the model has a wide variety of applications. Before considering 
Black-Scholes model, there is a number of assumptions that should be made. Fischer Black calls 
them “ideal condition” of the market (Black and Scholes (1973)). These assumptions are important 
to emphasize because it is well-known that stock markets are often volatile compared to other parts 
of the economy. 

There are five underlying assumptions: 

1. First assumption that should be made is information about values of short-term rates is available 
and short-term interest rates are constant (Black and Scholes (1973)). 

2. Secondly, stock pays no dividend (Black and Scholes (1973)). 

3. Thirdly, transaction costs that occur while buying or selling securities are eliminated (Black and 
Scholes (1973)). 

4. Fourthly, it is achievable to borrow fraction of price of stock one wants in order to hold the stock 
(Black and Scholes (1973)). 
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5. The last assumption states that short selling of a security in necessary situations is allowed in 
the market (Black and Scholes (1973)). 

Thanks to these assumptions, option price will be the function of the time period and stock price 
only. In the following paragraphs, option price will be reduced to the function of stock price only for 
simplicity. 

Generally, option values increase when stock prices rise. A positive relationship between option 
value and stock price may be easily seen from the following graph (Black and Scholes (1973)). As it 
can be seen from the figure, graphs representing the relationship between the option price and stock 
price at different time periods (T \, T 2 . T 3 ) lie below 45-degree line, which shows that option prices are 
more volatile than the stock prices (Black and Scholes (1973)). The volatility of option prices lead to 
the following statement: if the price of the stock increases by a certain amount, greater percentage 
change will be generated in option prices. The graph illustrated below shows what the paper seeks to 
explain through Black-Scholes model. 



Figure 1: The relation between option value and stock price 


2.2 Transformation into an Ordinary Differential Equation 

Black-Scholes equation is given by the following expression: 

dc i 2 2 dc 2 dc „ n 

-m + r s 8? + TS m- rC = 0 ' 

where C(s, t) = price of option, s = price of the stock, t = period of time, r = interest rate (Company 
et al. (2007)). Firstly, it is useful to transform this partial differential equation (PDE) into an ordi¬ 
nary differential equation (ODE) by proposing the following solution: C(s,t) = C(s)e xt . Given that 


dC 


= C(s) * \e xt and = 
dt K J ds 


dC dC{s) 


e xt , by substituting these equations into the PDE we get: 


ds 


C(s ) * Xe xt + — u 2 s 2 ^-j-^—^ * e xt + rs * e M _ r C(s ) * e xt = 0. 
2 as z ds 

The next step is to rearrange the equation to get second order ODE: 

'1 


„\t 




= 0 . 


The latter expression can be reduced to the following equation: 


^^) + „^ + C (.)( A _ r) _0 


since e 


A t 


7^0. 
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2.3 Euler equation 

To get rid of the coefficient of the first term lets divide everything by 1/2 ct 2 : 

2 d 2 C(s ) 2 r dC(s ) 2 (A — r) 


“ ds 2 + ^* S ~ ds 
This equation reminds us the Euler equation: 


+ 


< 7 - 


-C(s) = 0 . 


\ 2 d2 y , d y , o n 

L{y) = x ^ + ax ^ + py = 0- 

with real constants a and (3 (Boyce and DiPrima (2009)). In our case, a = % and j3 = 2 ^ r ) , which 
are positive constants. Euler equation has the solution of the form 

y = x ri + x r2 

in case of distinct real roots, and characteristic equation of the form: 

F(r) = r(r — 1) + ar + (3 = 0 


(Boyce and DiPrima (2009)). 


2.4 Solution of the Black-Scholes equation 

By the assumption given, a and r are positive real numbers because r is an interest rate and u 
is volatility of the stock as noted earlier in the paper. Now, a solution in the form of C(s) = s k can 
be proposed and applied to the Black-Scholes equation. The following derivations will be useful in 
solving our problem: 


C(s) = s\ 


ds 


d 2 C(s) 
ds 2 


= k(k — 1 ) * s 


k—2 


Substituting the derivations back into the earlier equation we get: 


-a 2 s 2 k(k — 1) * s k 2 + rsk * s k 1 + s k (A — r) = 0 . 


The next step is to take s k out of bracket and derive characteristic equation introduced earlier: 

= 0 . 


s k * 


^a 2 k(k — 1) + k(r — ^u 2 ) + (A - r) 


^cr 2 £; 2 + rk + (A — r) = 0. 


To find the roots of characteristic equation, let us find discriminant: 


D = (r — (jcr 2 ) 2 — 4 * ^<t 2 (A — r) = r 2 + ra 2 + — -2A<r 2 > 0 


T ' * ' 2“ v ' ' 7 . ' 4 

by assumption. So, the two distinct roots of characteristic equation will be: 


k\.2 = 


(r — \<J 2 ) ± \jr 2 + ra 2 + ^ — 2\a 2 


<X 2 . Q 

T * 2 


r 1 


r 2 + ra 2 + w- — 2Au 2 


a z 


Therefore, the solution of our problem can be written as: 

Jr 2 + ra 2 + ^--2\a 2 

C(s) = Cis^“ 2+ ^ + C 2 s — 2 " 


\Jr 2 +ra 2 + £- 2Xa 2 
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The solution represents option value as a function of stock prices. By the assumption ci and C 2 
must be positive constants because of a positive relationship between option price and the stock price 
introduced earlier. 


Simulations of Small, Medium, and Large-Cap Stock Prices 


Based on Parameters From Monthly Stock Returns 1926-2008 



Figure 2: Fluctuations in stock prices from 2000 to 2009 

The figure 2 above shows fluctuations in stock prices from 2000 to 2009 time period. The Black- 
Scholes presented in the paper is useful to explain, predict and estimate option prices based on stock 
prices in the financial world. Black-Scholes model gives more accurate estimates of option prices than 
other earlier developed models because it takes into account such factors influencing the stock prices as 
transaction costs, riskiness of assets, illiquid markets (Ankudinova and Ehrhardt (2008)). Therefore, 
the model is used to estimate European call options, which consolidates its role in applied economics 
(Barad (2014)). 

Black-Scholes model focuses on option or security that is held for a certain period of time and 
gives the owner right to make market operations such as buying and selling. Two types of securities 
can be specified: American options and European options (Black and Scholes (1973)). The difference 
of American options from European options is its quality that allows the owner to buy or sell the 
security until the maturity date, whereas the latter one does not allow conducting market operations 
until the security matures. 

According to empirical tests made by Black and Scholes, estimated option values deviate slightly 
from what they are in practice. Although those who demand stocks and bonds pay higher (to a 
small degree) prices for the products in securities markets, suppliers receive payments fairly close to 
what the formula calculates (Black and Scholes (1973)). This gap occurring between prices paid and 
received by demanders and suppliers may be understood by transaction costs — costs associated with 
exchanges — which occur as a result of a variety of services in an industry. 

3 The essence of the Black-Scholes Equation 

There are two types of options that can be specified: “American option” and “European option” 
(Black and Scholes (1973)). American options are the ones that are checkable on demand, particularly 
they can be returned before the maturity date. Whereas European options are the ones that can be 
returned only on a specific date, when they mature. So, American options are more liquid than 
European options. 

As noted earlier, an option is more valuable when the stock price is higher (Black and Scholes 
(1973)). Also, option depends on the maturity date, particularly the date of expiration. If maturity 
date is over a long period of time, then payments, particularly dividends that are paid on specific 
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periods of time for the option will be less. On the other hand, if maturity date is during a short period 
of time, then dividends are higher. 

Since Black-Scholes equation is a theoretical prediction of stock movements in the market, there 
are some restrictions that should be noted. Beside theory, there is a real world in which conditions 
in the stock market may not be that ideal as predicted by the model. Also, the reason theoretical 
models are built in finance and economics is that it is difficult to test or do an experiment with the 
real world. For example, in physics the experiments can be done inside the laboratories as well as in 
chemistry. However, experimenting with financial economics situations imposes high costs as they are 
closely connected with markets, which exist in our real world. So, making an experiment with the 
financial markets is difficult because of globality of problems. 

The only information for conducting policies or solving problems in finance is past data collection. 
Data in the past is similar with a tool that helps to understand general patterns of markets, movements 
of asset prices, behaviors of economic agents, relationships and co-movements of financial variables. 
With data collected in the past, it is possible to draw graphs relating the economic variables, showing 
their trends and movements. Also, data helps to predict future trend in markets, however only roughly. 

Now as we understand the importance of theoretical models, there are some conditions that should 
be specified regarding the Black-Scholes model. This is a complete list of assumptions additional to 
that noted earlier in the paper: 

1. There are perfect information about short - term interest rates and their movements are constant 
over long period of time (Black and Scholes (1973)). 

2. Movement of stock prices is random and its variance is proportional to the square of stock prices; 
therefore distribution of stock prices is lognormal over some finite period t (Black and Scholes 
(1973)). 

3. Unlike the real world, where stocks pay dividends to the shareholders, in the model stocks do 
not pay any payments (Black and Scholes (1973)). 

4. Only European options are to be considered in the model, so they are returned on a maturity 
date (Black and Scholes (1973)). 

5. Unlike in financial markets where transaction costs exist in all operations, it is assumed that 
buying or selling stocks do not impose any transaction costs (Black and Scholes (1973)). 

6 . Any fraction of the price of security can be borrowed at short — term interest rate (Black and 
Scholes (1973)). 

7. There is a possibility of short selling with no penalty or any costs. A seller will accept the price 
that buyer tells, agrees to meet at some specific time in the future, and pay the amount equal 
to that of the security price (Black and Scholes (1973)). 

Due to these assumptions, option value will depend only on time and stock price; and other 
variables are taken to be constant so as to simplify our model (Black and Scholes (1973)). Therefore, 
the value of the option w reduces to the following simple function: 

w = f(x,t ), 

where x is the price of the stock and t is the period of time (Black and Scholes (1973)). The expression 
above tells us how the value of the option changes if one of the variables, price of the stock or period 
of time, changes. The assumption that stocks pay no dividends gives an advantage of getting to more 
complicated problems with options. One example is under certain conditions, the formula can be 
applied to American options, which can be issued before maturity date (Black and Scholes (1973)). 
Black and Scholes (1973) studied this particular case in 1973. 

The Black-Scholes equation helps to calculate not only option value, but also more complicated 
assets such as warrants, which are liabilities of corporations (Black and Scholes (1973)). Warrants 
are generally considered as options. Also, the value of corporate liabilities may be calculated via the 
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formula. Often corporate liabilities are not viewed as options. We may consider a case of a company, 
which has assets of shares of another company (Black and Scholes (1973)). Also, bonds are “pure 
discount bonds”, which indicates a bond pays fixed amount of money (Black and Scholes (1973)). 
Assume maturity of 10 years. Also, assume that a company has a restriction of paying no dividend 
until the maturity date. Then, assume the company is planning to sell all of its stocks after 10 years 
and pay the amount for holders of bonds (Black and Scholes (1973)). So, these certain conditions let 
us to calculate a value of corporate liabilities using Black-Scholes equation as the assumptions make 
corporate liabilities similar to options. 

4 Partial Differential equation 

This part of our paper will greatly consist of the material from the book by Salsa (2008). Also 
note that Zhexembay and Pak (2016) focus on numerical solution of nonlinear Black-Scholes equation 
using Finite Element Method. Let us construct a differential equation that will help us to describe 
the evolution of V(s,t). The following hypotheses are established: 

• S follows a lognormal law; 

• The volatility a is constant and known; 

• There are no transaction costs or dividends; 

• It is possible to buy or sell any number of the underlying asset; 

• There is an interest rate r > 0, for a riskless investment. This means that 1-dollar in a bank at 
time t = 0 becomes e rT dollars at time T; 

• The market is arbitrage free. 

The last hypothesis is crucial in the construction of the model and means that there is no op¬ 
portunity for instantaneous risk-free profit. It could be considered as a sort of conservation law for 
money! We can translate this principle into mathematical terms through the notion of hedging and 
the existence of self-financing portfolios. The fundamental idea is to calculate the return of V through 
formula and then to build a riskless portfolio II. This portfolio contains shares of S and the option. 
II must increase at the current interest rate r, i.e. dll = rlldf, which turns out to coincide with 
the fundamental Black-Scholes equation. Now let us move to the calculation of the differential of V 
through the means of the formula. Since 


dS = fiSdt + adB , 


the obtained result is 

dV = [V t + aSV s + 1/2 ii 2 S 2 V SS \dt + aSV s dB. (1) 

In the formula (1) we have risk term aSVsdB. So, our next goal will be an elimination of this 
term. This can be acquired by establishing a portfolio II consisting of the option and a quantity —A 
of underlying: 

U = V- SA. 

This operation is valuable financial procedure called hedging. Now, let us turn out attention to the 
particular period of time, say (t, t + dt) during which II goes through a variation dll. If we succeed in 
keeping A equal to its value at t during the interval (t, t + dt), the variation of II is given by 

dU = dV - A dS. 

Since the mentioned formula is a cornerstone of the whole construction, it should be explained properly. 
Implementing 1 the we acquire: 

dU = dv - A dS = [V t + nSV s + ^a 2 S 2 V ss - nSA]dt + aS(V s - A )dB. (2) 
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Thus, if we choose 


A = V a , (3) 

where A is the value of V s at t, we eliminate the stochastic component in (2). The development of 
the portfolio II is now deterministic and its dynamics can be described by the equation: 

dn = [14 + l -a 2 S 2 V ss \dt. (4) 

The choice of (3) seems inexplicable, but it can be justified by the fact that V and S are dependent 
and the random component in their dynamics is proportional to S. Thus, we the linear combination 
of V and S is chosen wisely, such component should vanish. 

Now let us implement the no-arbitrage principle. Investing II at the riskless rate r, after a time 
dt we have an increment rlldt. Comparison between rlldt and dll is given by (4). If dll > rlldt, we 
borrow an amount II to invest in the portfolio. The return cfll would be greater of the cost rlidt, so 
that we make an instantaneous riskless profit 

dn — rlldt. 

if dn < rlldt, we sell the portfolio n investing it in a bank at the rate r. This time we would make 
an instantaneous risk free profit 

rndt — dn. 

Therefore, the arbitrage free hypothesis forces 

dn = [V t + T 2 S 2 V ss }dt = rlldt. (5) 

Substituting n = V — S A = V — V s S into (5), we obtain famous Black-Scholes equation: 

CV = Vt+ l -(J 2 S 2 Vss T rSVs -rV = 0. (6) 

Since the coefficient of V ss is positive, (6) is a backward equation. In order to get well-posed problem, 
we need to impose final condition (at t = T), a side condition at S' = 0 and one condition for S —> +oo. 

• Final conditions (t = T) 

Call. If at time T we have S > E then we exercise the option, with a profit S — E. If S < E, 
we do not exercise the option with no profit. The final payoff of the option is therefore 

C(S, T ) = max[S -E,0] = (S- E )+, S > 0. 

Put. If at time T we have S < E, we do not exercise the option, while we exercise the option if 
S < E. The final payoff of the option is therefore 

P(S, T) = ma x[E - S, 0] = (E - S)+, S > 0. 

• Boundary conditions (S = 0 and S —> +oo) 

Call. If S = 0 at a time t, S = 0 thereafter, and the option has no value; thus 

C(0,t) = 0,t > 0. 

As S —> Too, at time t, the option will be exercised and its value becomes practically equal to 
S minus the discounted exercise price, so 

C(S, t)-(S - e-^-^E) -T 0 as S -> Too. 

Put. If at a certain time is S = 0, so that S = o thereafter, the final profit is E. Thus, to 
determine P(0,t) we need to determine the present value of E at time T, that is 

P(S,t) = Ee~ r( ' T ~ t \ 

If S —> Too, we do not exercise the option, hence 

P(S, t) = 0 as S -> Too. 
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Solution of the Black-Scholes equation 

Let us summarize our model in the two cases. 


Black-Scholes equation 


1 


V t T =:<J 2 S 2 V ss T rSVs -rV = 0. 


(7) 


Final payoffs 


C(S,T) = (S — E) + 

(call) 


P(S,T) = ( E-S) + 

(put). 


Boundary conditions 



C(S,t) - (S - e-^-^E) -> 0 

as S —»• Too 

(call) 

P(0,t) = Ee~ r ( T "*), P(S,t) = 0 

as S —>• Too 

(put) 


The problem above can be simplified to a global Cauchy problem for the heat equation. Thus, the 
explicit formulas for the solutions can be obtained. Firstly, a change of variables should be performed 
so as to reduce the Black-Scholes equation to constant coefficients and to pass from backward to 
forward in time. Note also that 1 /cr 2 can be considered an intrinsic reference time while the exercise 
price E gives a characteristic order of magnitude for S and V. Thus, 1/cr 2 and E can be used as 
rescaling factors to introduce a dimensionless variable. Let us set 


, ,S 

x = In | — 


T = ^o- 2 (T-f), 


uj(x, t) = ^V(Ee x ,T — ^). 

L, 


E a z 

When S goes from 0 to = oo, x varies from — oo to +oo. When t = T we have r = 0. Moreover: 

Vt = —~a 2 Euj T 


_ E E E 

V s — —U] Xi Vs s — ~g2^ x ~g2^ xx ' 

If we put this into (7) we will end up with the following: 

1 2 1 

--cr u T + - {-u x T u xx ) + ru x - ru = 0 


or 


oj t = tjj xx + (k — l)u i x — ku 

where k = ^ is a parameter with no dimension. If we set 


k -1 (k+l) 2 

u(x,r) = e 2 x 4 T u(x, r), 


we find that v satisfies 

v T — u xx = 0, —oo < x < +oo, 0 < r < T. 

It is worth mentioning that the final condition for V is an initial condition for v. Performing 
following several steps we find out that 


v(x,0) = g[x ) 


for the call option, and 


v(x,0) = g{x) 


e lV+i)x _ e i(fe-i)x j x > o 
0 , x < 0 


e |(fe- l)x _ e \{k+ 1)® x < o 

0 , x > 0 
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for the put option. 

Now we can use the preceding results to derive the formula of the solution. The solution is unique 
and it is given by formula 


z '(x,t) = 


1 


g{y)e 


(x- v r 


dy. 


y/‘ 4vrr Jr 

To have a more general formula, let y = \/ 2 tz + x. Then, focusing on the call option: 


^ _ _ x 2 

v(x,t) = —j= / g(V2rz + x)e~^~ dy 

V 4ttt Jr 


V2^ 


/ oo _ poo 

e §( fe +l)(v / 2 tz+x)-±z 2 fa _ / e t(fc-l)(\/2TT+a;)-t2 2 ^ 

s lx 


After modifying those two integrals, we get 

v(x,t) = e^ k+1)x+ ^ k+1)2r N(d+) - e | ( k -i)x+U k -i) 2 T N (d_) 

where 

N(z) = —j= [ e~^ y dy 
V J —oo 

is the distribution of a standard normal random variable and 

d± = WT +l J k±1) ' / ^' 

Returning to the original variables, for the call we have: 

C(S,t) = SN{d + ) - Ee~ r ^ r ~ t ' ) N(d-) 


with 

The formula for the put is 
It can be shown that 


d± = 


In (S/E) + (r±^)(T-r) 


(J\jT — T 

P(S,t ) = Ee- r( ?- t) N(d-) - SN(d+). 

A = C s = N(d+) > 0 for the call 
A = P s = N(d+) — 1 < 0 for the put 

(Salsa (2008)). We should pay particular attention to the fact that both C s and P s are strictly 
increasing with respect to S. Thus, the functions C, P are strictly convex functions of S, for every t, 
namely C ss > 0 and P ss > 0. 

• Put-call parity. Put and call options with the same exercise price and expiry time can be 
connected by forming the following portfolio: 

U=S+P-C 

where the minus in front of C shows short position (negative holding). For this portfolio the final 
payoff is 

II(S, T) = S + (E - S) + - {S - E) + . 


If E > S, we have 
while if E < S 


n (5,T) = S + (E-S)-0 = E, 
11(5, T) = S + 0- {S - E) = E. 
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Therefore, at expiry, the payoff is always equal to E and it forms a riskless profit, whose value at t 
must be equal to the discounted value of E, since the no-arbitrage condition was imposed. So, we find 
the subsequent relation (put-call parity ) 

S + P - C = Ee-rV-Q. (8) 

Formula (8) reveals that, with the value of C (or P) available, the value of P (or C) can be obtained. 
From (8), since Ee^ r ^ T ~^ < E and P > 0, we get 

C(S, t) = S + P — Ee- r{T - t] > S — E 


and since C > 0, 


C(S,t) > ( S-E) + . 


It can be observed that the value of C is always greater than the final payoff. However, this property 
does not hold for a put. In fact, 


P(0,t) = £e“ r(r_f) < E 


so the value of P is less than the final payoff when S approaches 0, and it is greater just before 
expiring. The figures 3 and 4 demonstrate that. 



Figure 3: The value for the European call option 



Figure 4: The value function of the European put option 

• Different volatilities. The comparison between the value of two options with the different volatil¬ 
ities (7 1 and (i 2 can be conducted through the means of the maximum principle. Let us assume 
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that the exercise price and strike time are the same for the both cases, E being the exercise 
price and T— strike time. Another assumption that we allow is that g\ > a 2 . Denote the values 
of the respective call options C^, C^. With decreasing amount of risk the value of the option 
should decline as well. What we want to show is that 

C ( 1 ) >C (2) , 5 > 0, 0 < t < T. 

Let W = C « - . Then 

W t + l -a 2 S 2 W ss + rSW s - rW = ^(a 2 - a 2 )S 2 c£>. (9) 

with W(S,T) = 0, W(0,t) = 0 and W —> 0 as S — > +oo. It is obvious that (9) is a nonhomogeneous 
equation, with the right hand side being negative for S > 0, since > 0. We know that W is 
continuous in the half strip [0, +oo) x [0, T] and disappears at infinity, it reaches its global minimum 
at (So,to). 

We claim that the minimum is zero and cannot be obtained at a point in the intervals (0, +oo) X 
[ 0 , T). The equation is backward, so to = 0 is not considered. Assume that W(So,to) < 0 with 5o > 0 
and 0 < to < T. Thus, 

W t (So,t 0 ) = 0 

and 

W s (S o ,t o ) = 0, W ss (S 0 ,t 0 )>0. 

Substituting S = So, t = to into (9) we observe a contradiction. Thus, W = C ^ — C ^ > 0 for 
5 > 0, 0 < t < T. 

In 1972, empirical tests on call-options by Fischer Black and Myron Scholes were done (Black 
and Scholes (1973)). The results of the tests show that actual prices at which agents of the economy 
buy and sell options deviate systematically from the prices predicted by the Black- Scholes model 
(Black and Scholes (1973)). Options are bought at consecutively higher prices from those predicted 
by the model, whereas options are sold approximately at the prices predicted by the model (Black 
and Scholes (1973)). It should be noted that the difference in prices paid by option buyers are higher 
for lower-risk stocks than for higher-risk stocks (Black and Scholes (1973)). The latter point makes 
sense because low-risk stocks are always preferable than high-risk stocks since the probability that 
low-risk stocks generate large profits and do not default are higher. For option buyers, there are high 
transaction costs involved in the real world. This fact might explain why formula underestimates 
the price paid by option buyers because the model is built under assumption of no transaction costs. 
According to Black and Scholes (1973), getting into account the magnitude of the transaction costs in 
the market, misestimation of the prices does not indicate potential profit opportunities for speculators 
in the market. 


5 Separation of variables method 

This section considers the partial Black-Scholes equation of the form 

dC 1 2 2 d 2 C dC n 

¥b"8? +r T =0 ' 

The aim of the section is to introduce separation of variables method in order to solve the equation 
and find the general solution. Let the function C(s,t ) to be written in the following form: 

C(s,t) = S(s)T(t). 


Then, the following partial derivatives can be derived: 


dC 

~dt 


ST', 
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d 2 C 

dt 2 


= S"T. 


Substituting the derivatives back to the original equation, we get: 

1 


ST' + ± a 2 s 2 S"T + rsS'T - rST = 0. 


For simplicity, let = a and r = b for now. Dividing the equation by ST, 

2 S" , S' T' , 

as — + bs— + — —b — 0. 

Rearranging and equating to a constant (c > 0), we get pair of two ordinary differential equations: 

2 5" S' , T' 

as — + bs— = b — — = c. 


The first equation 


S S 


as 2 S" + bsS' = cS 


can be solved by Euler equation method introduced earlier. Dividing both sides of the equation by a 
and rearranging: 


Characteristic equation is therefore: 


or 


s 2 S” + b sS' --5 = 0. 
a a 


d(d — 1) + -d — - = 0 
a a 


d 2 + (-1 )d — 1 = 0. 


The characteristic equation has the following solutions (assuming D > 0): 


di,2 = 

Therefore, the general solution of the equation is 




, . (1-|) + a/(|-D 2 -4(-§) (1-|)-a/(|-1) 2 -4(-§) 

S(s) = s 2 + s 2 

Substituting respective expressions instead of a and b, we get the solution in explicit form: 


S{s) = s 2 +s 2 

Next, the aim is to find the solution of T(t) from: 

r 

b - — = c. 

T 

T - (6 - c)T = 0. 

So, the solution of the equation is therefore: 

T(t) = e (fe - c)t . 

Since the two solutions of the form S(s ) and T(.s’) are found, the final solution of the Black-Scholes 
equation 


dC 1 2 2 d 2 C dC 

-m + 2 r ‘ -W + T ’^ = 0 


can be expressed as 


C(s,t) = S(s)T(t ) = 




+ s" 


D (b—c)t 
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Figure 5: Call option 
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Figure 6: Put option 


• Constructing a plot 

Nowadays, there are lots of opportunities to see the graphical solution of the Black-Scholes equation. 
We decided to rely on one of them and show the plot of this equation. Statistic Online Computational 
Research (SOCR) allowed us to construct these two graphs. In the graph 5 the Exercise Price E was 
taken to be equal to 100, Interest Rate r = 0.5, Dividend Rate <5 = 0, Variance a = 0.3 and Time 
to Expiry T — t = 1. The graph was plotted with respect to the Stock Price S and Price V. On the 
graph 6 the values of the variables were assumed to be the same and the plot was made according to 
the Put option. 

6 Conclusion 

To conclude, Black-Scholes model is highly appreciated in quantitative finance because of its ac¬ 
curate and useful estimation of stock prices. Black-Scholes equation represents derivation of option 
pricing though taking into account such factors as time period t, risk-free interest rate r and volatility 
of stock prices a (Sheraza and Preda (2014)). Derived solution for the option value is closely related to 
corporate liabilities, therefore, the formula derived may be used to securities, including common stock 
and bond (Black and Scholes (1973)). This feature of Black-Scholes model illustrates its flexibility and 
efficiency of being applied to different contexts in the financial world. In this paper, we proposed a 
new method of solving the famous Black-Scholes Equation. Separation of variables method was used 
to derive a solution to the partial differential equation. 
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